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1 1 Many multiferroic materials, with various chemical compositions and crystal structures, 

have been discovered in the past years. Among these multiferroics, some perovskite 
manganites with ferroelectricity driven by magnetic orders are of particular interest. In 
^ these multiferroic perovskite manganites, not only their multiferroic properties are quite 

prominent, but also the involved physical mechanisms are very plenty and representative. 

I"""""'" In this Brief Review, we will introduce some recent theoretical and experimental progress 

on multiferroic manganites. 

o 
o 

1. Introduction 



Keywords: Multiferroicity; manganites; spiral spin order; exchange striction; CaMn70i2- 



Multiferroicity denotes the co-existence of more than one primary ferroic order pa- 
rameter simultaneously in a single material.^ These ferroic order parameters can be 
ferromagnetism, ferroelectricity, ferroelasticity, and ferrotoroidicity.^ To date, the 
most studied multiferroics are those materials with ferromagnetism (or antiferro- 
magnetism) and ferroelectricityP Since ferromagnetism and ferroelectricity are very 
useful in spintronics, information storages, sensors, etc., multiferroics with sponta- 
neous magnetization and polarization are promising candidates to design advanced 
devices with faster speeds, more functions, and energy saving. 

However, in many materials, the magnetism and ferroelectricity are mutually 
exclusive,^ which once blocked the development of multiferroicity in the last cen- 
tury. Till 2003, two milestone works, the discovery of magnetic-field-controllable 
ferroelectric polarization in TbMn03 crystals and a giant ferroelectric polarization 
in BiFe03 films, revived the research interests on multiferroicityEE Since then, 
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Fig. 1. (a) The perovskite structure. Center: B-site Mn cation; Corner: A-site rare earth/alkaline 
earth cations; The octahedral cage surrounding Mn is constructed by oxygen anions, (b) In 
many perovskite oxides, oxygen octahedra tilt and rotate collectively. These distortions reduce 
the nearest-neighbor hoppings (thus the bandwidth) of e g electrons, (a) and (b) are drawn by 
VESTAESl ( c ) Schematic drawing of split of Mn's 3d levels by the crystal field and Jahn- Teller 
distortions. 



multiferroics have become a flourishing research area. More and more multiferroic 
materials have been discovered and the understanding of underlying physical mech- 
anisms have also been pushed forward gradually. 

Up to date, there are already many experimentally verified multiferroic materi- 
als, most of which are oxides although there are also some exceptions, e.g. some flu- 
orides and organic tetrathiafulvalene-p-bromanilEHU According to Khomskii's clas- 
sification,^ some of them belong to the type-II multiferroics (or so-called magnetic 
multiferroics) , in which the origins of ferroelectric polarizations are relevant to some 
particular magnetic profiles, while in the type-I multiferroics the origins of polar- 
izations are almost independent of magnetism. Therefore, the type-II multiferroics 
have strong intrinsic magnetoelectric couplings, or more precisely, their ferroelectric 
polarizations can be significantly affected by tuning their magnetism, while in most 
type-I multiferroics, the magnetoelectric couplings are usually weak. Therefore, the 
type-II multiferroics are more interesting in physics, and will be also important in 
future applications. 

In this Brief Review, multiferroic perovskite manganites will be introduced, 
which are typical type-II multiferroic materials. Among various type-II multifer- 
roics, these perovskite manganites are the most important because not only their 
multiferroic properties are quite prominent, but also the involved physical mecha- 
nisms are very plenty. For example, the largest polarization of type-II multiferroics, 
always belongs to perovskite manganites so far, although this record has been up- 
dated year by year. In addition, almost all known physical mechanisms involved 
in the type-II multiferroicity, exist in perovskite manganites. In fact, most of these 
mechanisms were firstly revealed in perovskite manganites. 

In the past years, there were already many excellent reviews on the multifcr- 
roicityHHEflfH] For this reason, the present Brief Review will select some recent 
theoretical and experimental progress on multiferroic perovskite manganites, espe- 
cially those works after above reviews. Some key physical issues which have already 
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been discussed in above reviews will also be updated according to the latest dis- 
coveries. Some unsolved debates will also be presented for further discussion. 

2. Multiferroic phases in .RM11O3 

2.1. Origin of polarization in the spiral spin phase 

In the beginning, Kimura et al. noticed that the ferroelectric polarization in 
TbMn03 relates to the incommensurate-commensurate (lock-in) transition at 28 
K, below which the magnetic modulation wave vector of Mn is locked at a constant 
valueP Then, by using the neutron diffraction technique, Kenzelmann et al. dis- 
covered that the Mn's magnetic structure of TbMn03 below 28 K is a transverse 
incommensurate spiral order, which breaks patial inversion symmetry and induces 
the magnetoelectricityPS Soon after, Arima et al. confirmed the spiral order in 
Tbi-zDy^MnOa, which's spiral period shrinks with x!^ 

Theoretically, in 2005, Katsura, Nagaosa, and Balatsky (KNB) first proposed 
a theory based on the spin supercurrent to explain the magnetoelectric effect in 
noncollinear magnets.^ Accroding to the KNB theory, a noncollinear spin pair 
gives rise to a ferroelectric polarization in proportion to ejy x (Si x Sj), where &ij 
is the unit vector pointing from site i to site j, and Si denotes the spin at site i. 
Later, Sergienko and Dagotto proposed that the (inverse) Dzyaloshinskii-Moriya 
interaction is responsible for the polarization in i?Mn03 (R denotes a rare earth 
cation, here is Tb or Dy)P^ The formula of Dzyaloshinskii-Moriya interaction is 
Dij ■ (Si x Sj) and is perpendicular to e^EUI 2 ^ Besides, Mostovoy derived a 
similar expression for polarization induced by incommensurate spin density wave 
states (in fact the cycloid spiral spin states) from the phenomenological Ginzburg- 
Landau approach.^ The KNB theory and Dzyaloshinskii-Moriya scenario seem to 
be equivalent formally, both of which depend on the spin-orbit coupling and are in 
proportion to the cross production of noncollinear spin pair (Si xSj). These theories 
were further confirmed by neutron diffraction measurements on Gdo.yTbo.sMnO;? 
which established the correlation between the spin-helicity and electric polariza- 
tion.^ 

However, in the KNB theory the polarization is purely electronic, while it is from 
ionic displacements in the scenario of Dzyaloshinskii-Moriya interaction. Since the 
(saturated) polarizations of i?MnC>3 (i?=Tb, Dy, ~E\ii_ x Y x ) are very weak (typically 
in the order of 0.1 /xC/cm 2 , which is already quite prominent among magnetic mul- 
tiferroics) comparing with traditional good ferroelectric materials, e.g. about 0.1% 
of BiFcC>3 (in the order of 10 2 /j,C/cm 2 ), it is not easy to direct measure the am- 
plitude of ionic displacements in i?MnC>3. Therefore, ab initio calculation becomes 
an alternative approach to judge which one (electronic or ionic) is the dominant 
contribution in i?MnC>3. The density functional theory calculation on TbMnC>3 by 
Malashevich and Vanderbilt found that the polarization from electronic contribu- 
tion was very weak (32 ^C/m 2 , far below the real value ~ 600 /zC/m 2 ) and with 
an opposite sign against the KNB equationPSl Only when the ionic displacements 
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Fig. 2. Schematic drawing of the magnetic moments in Tt>o.4iDyo.59Mn03 projected onto the be 
plane at (a) 15 and (b) 30 K. The trajectories of the Mn spins are also shown, (c) Temperature 
dependence of lattice modulation wave number (q^, which is twofold of the magnetic modulation 
wave number QMn) in some Tbi-^Dy^MnOs. The ferroelectric Curie temperatures are indicated 
by arrows. Repri nted figure with permission granted from |T. Arima et al, Phys. Rev. Lett. 96| 
(2006) 097202 B3 Copyright ©(2006) by the American Physical Society. 

were taken into account, a realistic polarization was reached in their calculation. 
Similar result was also obtained in Xiang et al.'s density functional theory workPS 
In short, these theoretical calculations supported the ionic displacement scenario, 
at least for TbMnC>3. However, in a following work, Malashevich and Vanderbilt 
changed their tune a little because they found the proportions of electronic/ionic 
contribution seriously depended on the spiral plane (ab or be) as well as the oc- 
tahedral rotations!^ In 2011, by using a diffraction technique which exploited the 
interference between the non-resonant magnetic scattering and the charge scattering 
arising from the ionic displacements, Walker et al. measured the ionic displacements 
in TbMnOa, which obtained very high precision values (up to several femto-meter) 
and decisively supported the microscopic models that attribute polarization to ionic 
displacements.^ 

2.2. Origin of the spiral spin order 

Above theoretical and experimental works have revealed how spiral magnetic orders 
generate ferroelectric polarizations. Then a following question is how to generate 
a spiral magnetic order? Or more specificly, why are there spiral spin orders in 
TbMnOa and DyMn03, instead of in LaMnC>3 or H0M11O3. This question is quite 
nontrivial since the underlying physical mechanisms can guide us not only to tune 
the multiferroicity in manganites and related materials but also to search for new 
multiferroics with better performance. 

In doped manganites, there are many phases and the competition between these 
phases may give rise to dramatic responses to external stimulations, e.g. the colos- 
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Fig. 3. Experimental magnetic phase diagram of BMnC>3. The horizontal axis stands for the de- 
creasing Mn-O-Mn bond angle (increasing oxygen octahcdra distortions) or namely the decreasing 
size of R. The Neel and lock-in transition temperatures are marked as open and closed triangles, 
respectively. Upper inset: Temperature dependencies of the wave numbers of lattice (<5;) or mag- 
netic (<5 m ) modulation. The A-type and E-type antiferromagnetic structures are shown in their 
regions. The gray area is the spiral spin order state. Reprinted figure with permission granted 
from |T. Goto et al, Phys. Rev. Lett. 92 (2004) 257201^ 1 Copyright ©(2004) by the American 
Physical Society. 



sal magnetoresistive effect and colossal electroresistive effect In fact, even in 
undoped i?Mn03, there are also many competing phases. When R is large (e.g. La 
and Pr) the magnetic ground state is A-type antiferromagnetic phase. 38 When R is 
very small (e.g. Ho and Tm), it becomes E-type antiferromagnetic order. 39 In the 
middle region (e.g. Tb and Dy) between the A-type and E-type antiferromagnets, 
the ground state is the spiral spin order, as shown in the Fig. [3| 32 

The first theoretical investigation was done as early as the discovery of multi- 
ferroicity in TbMn03. In 2003, Kimura et al. studied a two dimensional frustrated 
classical spin Heisenberg model with ferromagnetic nearest-neighbor exchange (Ji), 
ferromagnetic/antiferromagnetic next-nearest- neighbor exchange along the a/b axis 
(J3/J2) respectively.^ The mean field phase diagram is shown in Fig. [4J With in- 
creasing |J 2 |, the magnetic transition occurs from the A-type antiferromagnet to 
the E-type antiferromagnet through the incommensurate structure. 

However, this classical spin model needs a strong antiferromagnetic next- 
nearest-neighbor exchange J2 to stabilize the spiral spin order and E-type anti- 
ferromagnetic order, like Ji = —J\ to obtain the q — 1/6 spiral. This abnormal 
large J2 seems to be unreasonable in the orthorhombic lattice because the ex- 
change path between next-nearest-neighbor sites is more complex and longer than 
the nearest-neighbor one. Furthermore, in this phase diagram, the incommensurate 
spin modulation fades away when temperature approaches zero, which disagrees 
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Fig. 4. Finite temperature (vertical axis) phase diagram of the two dimensional J± — J2 — J3 
model. Here J3/J1 = 0.01. q is the wave number of spin structure along the qseudocubic x and 
y axes, q = denotes the A-type antiferromagnetic order (ferromagnetic in the two dimensional 
lattice) and q = 1/4 denotes the E-type states. The middle q's stand for the incommensurate 
spin modulation. At that time, the spiral spin order had not been revealed. Reprinted figure 
with permission granted from |T. Kimura et al, Phys. Rev. B 68 (2003) 060403(R) 02l Copyright 
©(2003) by the American Physical Society. 



with the experimental phase diagram. 

Mochizuki and Furukawa extended the frustrated classical spin Heisenberg 
model to the three dimensional lattice, including more interactions such as the 
single-ion anisotropy, Dzyaloshinskii-Moriya interaction, and cubic anisotropyESE?] 
With proper parameters, the experimental phase diagram can be well reproduced, 
as shown in Fig. [5] The main conclusion is that the next-nearest-neighbor spin 
exchanges enhanced by the orthorhombic lattice distortion (the so-called GdFeC>3 
type distortion) induces the spiral spin order, while the spiral plane (ab vs. be) is 
exquisitely controlled by tuning the competition between the single-ion anisotropy 
and Dzyaloshinskii-Moriya interaction. 

Using this model, Mochizuki et al. studied many physical properties of multifer- 
roic i?MnC>3, such as the reorientation of polarizationJHllillelectromagnonsJ^ mag- 
netostrictions,^ picosecond optical switching of spin chiralityj^S magnetic switching 
of ferroelectricityPSl For more details of Mochizuki et al.'s model, readers can refer 
to their recent long paperPS 

Besides these models based on pure classical spins, there are also quantum 
models for i?MnC>3, e.g. the double-exchange model EHIMtlll In the double-exchange 
model, e g electrons are treated as itinerant fermions, while t 2g electrons are localized 
which's spin texture provides a magnetic background.^Hm The double-exchange 
model has been extensively studied in the past decade to understand phase compe- 
titions and the colossal magnetoresistive effect in manganites bulks, 55 56 as well as 
in thin films and heterostructuresf^22l which has been proved to be quite successful 
for manganites 12211^1 

For undoped i?MnC>3, both two e g orbitals (d x 2_ y 2 and d 3z 2_ r 2) must be con- 
sidered in the double-exchange process. Beside the ferromagnetic double-exchange, 
there are antiferromagnetic superexchanges between neighbor t% s spins. In addition, 
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Fig. 5. A typical finite temperature magnetoelectric phase diagram of the Mochizuki-Furukawa 
model as a function of J2. AFM + WFM is the A- type antiferromagnetic phase with weak ferro- 
magnetism due to the spin canting. PE/FE denotes paraelectric/ferroelectric phase, respectively. 
Reprinted figure with permission granted from M. Mochizuki et al., Phys. Rev. B 80 (2009) 
|134416| ^I Copyright ©(2009) by the American Physical Society. 



the Jahn- Teller distortions also have to be taken into account. 

In an early study, Hotta et al. investigated the phase diagram of i?Mn03 using 
this two-orbital double-exchange model!^ However, only the A-type and E-type 
antiferromagnetic phases were reproduced, without the spiral one, as shown in 
Fig.[6]ja). Then, in Sergienko and Dagotto's work, the spiral spin order was obtained 
when the Dzyaloshinskii-Moriya interaction was included, as shown in Fig.[6jb)p2l 
However, the required Dzyaloshinskii-Moriya interaction is too large, far beyond 
the real intensity. 

To overcome these problems, Dong et al. included the next-nearest-neighbor 
superexchange J 2 into the double-exchange model (without the Dzyaloshinskii- 
Moriya interaction since it is very weak)Psl Due to the GdFeC>3 type distortion 
of orthorhombic lattice, the next-nearest-neighbor exchange along the 6-axis (in 
the Pbnm notation) J 2 b is stronger than J 2a along the a-axis. Dong et al.'s calcula- 
tion found that a weak J 2 b was enough to stabilize the spiral spin order. However, 
to obtain the realistic spiral periods, e.g. ~ 7.2 cells for TbMnC>3 and ~ 5.2 cells 
for DyMnC>3 (along the qseudocubic x and y axes), the contribution from Jahn- 
Teller distortions is essential. In these orthorhombic perovskites, the lattice constant 
along the c-axis shrinks, which corresponds to the Jahn- Teller Q 3 mode, and the in- 
plane (ab plane) distortion (Jahn- Teller Q 2 mode) is also prominent. For undoped 
i?MnC>3, these is an approximate ratio \Q 2 \ ~ — V^Qs, which induces the special 
d 3x 2_ r 2 /d 3y 2_ r 2 type of orbital orderf^EH With proper Jahn- Teller distortions, a 
tiny J 2 b, in the order of 10% of the nearest-neighbor superexchange Jaf, is enough 
to obtain the realistic spirals, as shown in Fig. [7^a, b). Such a weak J 2 b and even 
weaker J 2a seems much reasonable in the orthorhombic lattice. 

With increasing Jaf and J 2 b, the ground state changes from the A-type anti- 
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Fig. 6. Ground-state phase diagrams of the two-orbital double-exchange model for _RMnC>3. The 
horizontal axis is the nearest-neighbor supercxchange. (a) A mean-field result on a 4 X 4 X 4 lattice. 
The vertical axis is the Jahn- Teller coefficient. FM denotes ferromagnetic order. OO (OD) denotes 
orbital order (disorder). A, C, E, G are various antiferromagnetic orders. Reprinted figure with 
permission granted from |T. Hotta et at, Phys. Rev. Lett. 90 (2003) 2472031 311 Copyright ©(2003) 
by the American Physical Society, (b) With the Dzyaloshinskii-Moriya interaction (vertical axis). 
Dashed lines are the low temperature Monte Carlo results on a 8 X 8 lattice. Solid lines are 
calculated in the thermodynamic limit. Here the ICM-FE denotes the (incommensurate) spiral 
order phase with a ferroelectric polarization. Reprinted figure with permission granted from I. A. 
|Sergienko et at, Phys. Rev. B 73 (2006) 094434^ 1 Copyright ©(2006) by the American Physical 
Society. 



ferromagnet to spiral phase, and finally to the E-type antiferromagnet, as found 
in real i?Mn03 with deceasing size of R. The following calculation on the three 
dimensional lattice also confirmed this A-S-E phase transition!^ 

The above zero-temperature phase diagrams should be further checked at finite 
temperatures. Monte Carlo simulations were employed to verify these spin orders 
and their phase transitions.^ The results are shown in Fig. [8} At low temperatures, 
all characteristic peaks of spin structure factor are prominent, suggesting robust 
spin orders. And the finite-temperature phase diagram qualitatively agrees with 
the experimental one. 

Later, Kumar et al. argued that even the weak J 2 was not necessary when using 
a finite Hund coupling (between the t 2g spins and e g spin) Their phase diagram 
was obtained by comparing zero-temperature energies of some candidate phases, 
which have not been further confirmed by unbiased Monte Carlo simulations. Other 
unexpected spin orders may emerge somewhere. In fact, the role of a weak J 2 b is not 
only to stabilize the non-collinear spin pairs, but also to make these cantings form 
uniform long-range chirality. Therefore, without this weak J 2 b, a uniform spiral may 
be not stable against canting states with random local chirality. 
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Fig. 7. Ground-state phase diagrams of two dimensional two-orbital double-exchange model for 
i?Mn03 with the next-nearest-neighbor superexchange Ji and Jahn- Teller distortions. Vertical 
axis: the nearest-neighbor superexchange; Horizontal axis: the next-nearest-neighbor superex- 
change along the fe-axis. A, S, E: A-type, spiral, E-type antiferromagnets, respectively, (a) J^a = 0. 
(b) J2a = 0.5J26- (c) The corresponding wave number (along the qseudocubic x and y axes) pro- 
file of magnetic phases in (b). (a-c) Re printed figure with permission granted from |S. Dong e1\ 
\al, Phys. Rev. B 78 (2008) 155121^ 21 Copyright ©(2008) by the American Physical Society, 
(d) The density of states for the A-type, spiral, E-type antiferromagnetic phases. Reprinted fig- 
ure with permission granted from S. Dong et al., Eur. Phys. J. B 71 (2009) 339^21 Copyright 
©(2009) by EDP Sciences, Societa Italiana di Fisica, Springer- Verlag. Here the energy unit is the 
double-exchange hopping amplitude (~ 0.2 eV). 



2.3. Multiferroicity of the E-type antiferromagnet 

In .RM11O3 , when the size of rare earth cation R is further reduced beyond Dy, the 
orthorhombic structure becomes unstable against the hexagonal structure. Even 
though, metastable orthorhombic structure remains available by using some special 
synthesis methods, e.g. under high-pressure conditions. In early days, experimental 
physicists were only interested in those spiral spin manganites like TbMnOa and 
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Fig. 8. (a) Typical spin structure factor of the four Monte Carlo simulations (corresponding to the 
four sets of parameters (asterisks) in Fig.JTJa)) at a low temperature, (b) Temperature-dependent 
polarization (P) and long-range spin correlation (C(L/2, L/2)) corresponding to the spiral phase 
with a wave number q = 1/6. All values are normalized to their saturation values, (c) Sketch of 
the finite-temperature phase diagram according to the transitions obtained in the Monte Carlo 
simulations. Reprinted figure with permission granted from |S. Dong et al, Phys. Rev. B 78 (2008)| 
|155121| 421 Copyright ©(2008) by the American Physical Society. 



DyMn03, while the E-type antiferromagnets (e.g. H0M11O3) were rarely concerned. 

In 2006, Sergienko, §en, and Dagotto predicted a large ferroelectric polarization 
driven by the spin zigzag chain in E-type antiferromagnets.^ In these systems, the 
mechanism for ferroelectricity is the exchange striction effect (relevant to Si ■ Sj), 
as shown in Fig. |9j which is independent of the weak spin-orbit coupling. There- 
fore, the induced polarization in orthorhombic HoMnOa can be greatly enhanced 
comparing with those spiral magnetic multiferroics. The theoretical value of polar- 
ization was estimated to be ~ 2 /iC/cm 2 in orthorhombic HoMnOa, 30 times more 
than TbMn0 3 P 

The following density functional theory calculations further confirmed the giant 
ferroelectric polarization (up to 6 /iC/cm 2 ) in HoMn0 3 The giant ferroelectric 
polarization in HoMn0 3 shows dual nature: a large portion arises due to quantum- 
mechanical effects of electron orbital polarization, in addition to the conventional 
polar atomic displacements.^ 

However, the first measurement of polarizations of HoMn03 and its analogous 
YMn0 3 only found moderate values: ~ 90 fiC/m 2 for HoMnC>3 and ~ 250 fiC/m 2 
for YMn03, which were far below the theoretical predictions.^ Even considering 
the polycrystalline factor, the expected values of corresponding single crystalline 
samples were not very exciting. 

The divergence between the theoretical predictions and experimental measure- 
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Fig. 9. Ferroelectric polarization in HoMnOs. (a) The original (nonpolar) configuration of a Mn- 
O-Mn bond, ipo is the bond angle, (b) A Monte Carlo snapshot of the E-type antiferromagnetic 
phase. The ferromagnetic zigzag chain links are shown as solid lines. The displacements of the 
oxygen atoms are also exaggerated, (c) The local arrangement of the Mn-O bonds. Left: without 
any magnetic order; Right: with the E-type antiferromagnetic order, (a)-(c) Reprinted figure with 
permission granted from |l. A. Sergienko et al., Phys. Rev. Lett. 97 (2006) 227204[ M1 Copyright 
©(2006) by the American Physical Society, (d) A typical charge density isosurface plot of the 
Mn02 plane. Arrows denote the direction of spins. Reprinted figure with permission granted 
from |S. Picozzi et al, Phys. Rev. Lett. 99 (2007) 227201^ 1 Copyright ©(2007) by the American 
Physical Society. 



merits has been gradually solved by improving samples' quality. Ishiwata et al. syn- 
thesized a series of high quality polycrystalline samples of i?Mn03 (i?=Eui_ x Y a; , 
Yi-yLuy, Dy, Ho, Er, Tm, and Yb), which showed a polarization up to ~ 800 
/xC/m 2 for those E-type antiferromagnetsF^ By multiplying a calibration factor 
(= 6), the value for corresponding single crystals was estimated to be ~ 5000 
/iC/m 2 , as shown in Fig. 10 a), which was much higher than the early results al- 
though still lower than the theoretical predictions. A possible reason is that the 
electric coercive forces of E-type antiferromagnetic multiferroics are very large, 
which require higher poling electric fields to obtain the saturated polarizations. By 
using a larger poling field (37.5 kV/cm -1 ), Pomjakushin et al. obtained a larger 
polarization up to 1500 ^C/m 2 in polycrystalline TmMn03, which seemed to be 
unsaturated yet P3 Recently, Ishiwata et al. grew YMn03 single crystals under high 
pressure.^ The polarization with 10 kV/cm _1 poling field reached 2200 /xC/m 2 at 
2 K, which did not saturate either.^ 

Comparing with bulks, thin films are more suitable to test higher poling fields. 
Nakamura et al. fabricated a monodomain single-crystal film of orthorhombic 
YMn0 3 on the (0 1 0) YA10 3 substrateP The measured in-plane polarization 
(with a poling field 10 kV/cm -1 ) reaches 8000 zxC/m 2 at 4 K, which seems to be 
close to the saturation. The following resonant soft x-ray and hard x-ray diffrac- 
tion revealed that the coexistence of E-type antiferromagnetic and cycloidal states 
below 35 KP The large polarization mainly comes from the E-type component 
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in YM11O3. Since this YMnC>3 thin film is not a pure E-type antiferromagnet, it 
is hopeful to expect a higher polarization in prototype E-type antiferromagnets at 
low temperatures. 

Also recently, Lee et al. synthesized rodlike large single crystals of orthorhom- 
bic HoMnC>3 utilizing the conventional Bi2C>3 flux method.^ Unexpectedly, their 
HoMn03 showed an incommensurate magnetic modulation (fc ~ 0.4) at all temper- 
atures,^ while an ideal E-type antiferromagnetic order should have a commensu- 
rate magnetic modulation (fc ~ 0.5) as found in previous polycrystalline samples.^ 
As a result of this incommensurate magnetic modulation, the exchange striction 
within the Mn's zigzag chains does not contribute a net polarization any more. 
However, another exchange striction which is between the Mn's and Ho's spins, 
contributes a ferroelectric polarization along the c-axis, which reaches 1500 /j,C/m 2 
at 2 K.E21 Currently, it is unclear why their single crystals are so different from pre- 
vious polycrystalline specimens. Lee et al. proposed that the presence of defects or 
residual strains might play an important role. Further experiments on orthorhom- 
bic HoMnOs single crystals are needed to solve above puzzles and verify the real 
saturated polarization!^! 

2.4. Ferroelectricity driven by dual magnetic effects 

In early studies, the Dzyaloshinskii-Moriya scenario (or the KNB theory) was 
applied to TbMnOs and DyMnC>3, while the exchange striction was applied to 
HoMnOs and YMnC>3. However, as shown in the above subsection, the latest ex- 
periments on YMn03 films and HoMnC>3 crystals showed some unusual results, 
which challenged the traditional viewpoint EHIZH 

According to Dong et aVs phase diagram, the phase transition between the 
spiral spin order and E-type antiferromagnetic order is first-orderPsl According 
to the experience in doped manganites with the colossal magnetoresistive effect, 
phase separations often emerge around first-order transition boundaries when some 
perturbations (e.g. quenching disorder) present pIHUl 

Therefore, it is reasonable to expect the phase separation between the spiral 
spin order and E-type antiferromagnetic order in i?Mn03 by fine tuning some con- 
ditions, e.g. the average R size. If so, the ferroelectricity can have dual origins: 
the Dzyaloshinskii-Moriya interaction effect as well as the exchange striction effect. 
The coupling between these two mechanisms may combine the giant ferroelectric 
polarization of exchange striction effect and the sensitive magnetoelectric response 
of Dzyaloshinskii-Moriya interaction effect in a single material. 

In 2009, Lu et al. synthesized Tbi-^Ho^MnOs polycrystalline samples, which's 
multiferroic transition temperatures showed a V-shape behavior with increasing 
X El] This V-shape behavior is a common symptom for phase separation tendency in 
the colossal magnetoresistive manganites.^ Meanwhile, Ishiwata et al. synthesized 
i?Mn03 (R=Eui- x Y x and Yi-^Lu,,) polycrystalline samples, and also found the 
V-shape behavior between the spiral spin phase and E-type antiferromagnets, as 
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Fig. 10. (a) Contour plot of polarization in the phase diagram of orthorhombic i?Mn03 with 
nonmagnetic R (Eui_ I Y I and Yi-^Luy). The polarizations in the E-type antiferromagnetic 
region are the polycrystalline values multiplying a calibration factor (= 6). Magnetic-field de- 
pendence of normalized polarization at 2 K for (b) Euo.6Yo.4Mn03, (c) Euo.25Yo.7sMn03, (d) 
Euo.iYo.gMnOs, and (e) LuMnOs. Reprinted figure with permission granted from S. Ishiwata et 
\al, Phys. Rev. B 81 (2010) 10041 1(R)| SSI Copyright ©(2010) by the American Physical Society. 



shown in Fig. |10[ a). The magnetoelectric responses (sensitivity) are significantly 
enhanced in the phase boundary region: Euo.25Yo.75Mn03 (located at the bottom 
of V-shape valley in Fig. [To][a)), as compared in Fig. 10 'b-e). This enhancement of 
magnetoelectric responses may arise from the phase separation between the spiral 
magnet and E-type antiferromagnet. As aforementioned, recent resonant soft x-ray 
and hard x-ray diffraction confirmed the coexistence of E-type antiferromagnetic 
and cycloidal states below 35 K in YM11O3 thin films.^1 

Despite the coexistence of E-type antiferromagnetic and cycloidal states, the 
presence of A-site magnetism may also provide one more source of ferroelectric 
polarizations, e.g. in the forementioned HoMnOa single crystals.^ The exchange 
coupling between the A-site rare earth and B-site Mn can generate a polarization 
via the exchange striction effect. In fact, even in DyMnC>3, which was previously be- 
lieved to be a prototype cycloidal spin multiferroic material, the exchange between 
Dy 3+ and Mn 3+ is responsible to the suppression of polarization at low temper- 
atures. Furthermore, recent experiments on Dyi-^Y^MnOs and Dyi_j,Ho a Mn03 
by Zhang et al. revealed that the a considerable portion of polarization in DyMn03 
is from the exchange striction effect between Dy 3+ and Mn 3+ P3EH Then it be- 
comes understandable why the improper polarization in DyMnC>3 (~ 2500 /xC/m 2 
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Fig. 11. Cartoon for the dual origins of ferroelectric polarization in Dyi-yHo^MnOa. The magni- 
tude and orientation of these spin arrows are only for guide of eyes. Upper: pure DyMnOa; Lower: 
Dyi—yHoyMnOs. Left: within the temperature region pb y , Tpp], where Tpp is the transition 
temperature of Mn's spiral spin order and Tp> y is the transition temperature of Dy's independent 
spin order. In this temperature region, Dy's m agnetic moments form an incommensurate modula- 
tion following Mn's modulation periodE3Z3 Right: at low temperatures (< 7b y ). Both the Mn's 
spiral order and Mn-Dy(Ho)'s exchange striction contribute to the polarization in cases (a), (c), 
and (d). In case (b), only the Mn's spiral contributes the polarization since the independent Dy 3+ 
magnetic order cancels the exchange striction between Mn 3 + and Dy 3 +. Reprinted figure with 
permission granted from |N. Zhang et al, Appl. Phys. Lett. 99 (2011) 102509[ ^ Copyright (2011) 
©by the American Institute of Physics. 
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Fig. 12. (a) Jt (the next-nearest-neighbor exchange) dependence of ground state polarizations. 
Ps comes from the exchange striction effect while Pas comes from the Dzyaloshinskii-Moriya 
interaction. Both Pg and Pas contribute to the polarization in the ab-spiral and E-type anti- 
ferromagnetic region. Reproduced from Ref. Reprinted figure with permission granted from M. 
|Mochizuki et al, Phys. Rev. B 84 (2011) 144409[ ^1 Copyright ©(2011) by the American Physical 
Society. 



under magnetic fields is much larger than other spiral magnets (typically be- 
low 1000 ^(C/m 2 )P^ The exchange between R and Mn also gives rise to prominent 
magnetoelectric responses in those E-type antiferromagnets with magnetic R, in 
contrast to the weak magnetoelectric responses in those E-type antiferromagnets 
with nonmagnetic i?P^ 
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In addition, there is one more case which can host the dual magnetic multifer- 
roic effects. Mochizuki et al. found the deformation of spin spiral can enable the 
exchange striction effect in those a6-plane spiral spin manganites, while in the fee- 
plane spiral spin manganites, this effect is canceled between different ab layersEHEH 
Mochizuki et al. also argued that there were noncollinear deformation of spins in 
the E-type antiferromagnets which was mostly believed to be collinear. If so, the 
Dzyaloshinskii-Moriya interaction would also contribute a small portion of polar- 
ization in the E-type antiferromagnets H1H1ZI However, this noncollinear deformation 
of spins needs further experimental verification. 

3. New multiferroic phases in doped manganites 
3.1. Multiferroic SOS phase 

Despite their strong magnetoelectric coupling, the main weaknesses of i?Mn03 are 
the low multiferroic temperatures (most below 40 K) and weak ferroelectric polar- 
izations. Therefore, one of the most important issues for magnetic multiferroics is 
to improve the multiferroic performance, including both the working temperatures 
and polarizations. 

In 2008, Kimura et al. reported that CuO is a "high" -temperature magnetic 
multiferroic material, which's multiferroic working temperature is 213 — 230 K.^In 
2010, Kitagawa et al. reported that Z-type Sr3Co2Fe2404i shows a spiral magnetic 
ordered state and a resultant magnetoelectric effect at room temperature.^! In 2011, 
Lee et al. reported that SmFeCVs magnetic multiferroic transition temperature 
reaches ~ 670 KP^ However, in all these new multiferroics, the polarizations remain 
quite weak (typically in the order of 100 /zC/m 2 for single crystals). 

Meanwhile, researchers are also searching in manganites for magnetic multifer- 
roics with better performance. In 2009, by using density functional theory, Gio- 
vannetti et al. predicted that half-doped manganites like Lao.sCao.sMnOs could be 
multiferroic materials with very prominent polarizations up to several /zC/cm 2 !^! 
This prediction based on an early idea of the site-centered vs. bond-centered charge 
order in half-doped manganites!^ However, till now, there is no direct experimental 
evidence to confirm the multifcrroicity in half-doped manganites. 

Also in 2009, Dong et al. re-investigated the phase diagram of quarter-doped 
two-orbital double-exchange model for manganites Although this phase diagram 
had once been studied by Hotta et aZjUE a new phase called SOS (for "spin 
orthogonal stripe" as explained below) was discovered by Dong et al. in the narrow 
bandwidth but weak Jahn- Teller distortions region, as shown in Fig. |13| In the 
narrow bandwidth and strong Jahn- Teller distortions region, it is the C!/4E 3 /4 
phase predicted by Hotta et aZEEHThe spin patterns of SOS and Cx^E^^ phases 
are shown in Fig. [I4^a,b). 

The C 1( /4E 3 /4 phase consists of zigzag chains of collinear spins, which can be 
viewed as a superstructure of 1/4 C-type antiferromagnet and 3/4 E-type antifer- 
romagnet. It belongs to a general family as C x Ei_a;l311llll The most studied case of 
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Fig. 13. Ground state phase diagram of quarter-doped two-orbital double-exchange model on the 
two-dimensional lattice. SOS is multiferroic. The horizontal axis Jaf is the nearest-neighbor su- 
perexchange and the vertical axis A is the Jahn- Teller coefficient. Full dots are those couplings 
where variational results have been confirmed by both the low-temperature Monte Carlo simu- 
lation and the zero-temperature relaxation on a 8 X 8 cluster. Open dots are cases where the 
Monte Carlo simulation does not provide a clear answer due to strong competition of metastable 
states. Reprinted figure with permission granted from |S. Dong et al, Phys. Rev. Lett . 103 (2009) 
I|l07204pl] Copyright ©(2009) by the American Physical Society. 

C x Fji- x is the CE phase (C1/2E1/2) in the half-doped manganites 

Different from the C1/4E3/4 phase, the SOS order consists of stripe blocks. 
Within each stripe, spins are collinear, which form the E-type antiferromagnetic 
order. But between nearest-neighbor stripes, spins are orthogonal to each other. 
This noncollinear spin pattern is not caused by neither the next-nearest-neighbor 
exchange frustration nor the Dzyaloshinskii-Moriya interaction. Instead, it origi- 
nates from the electronic self-construction at the particular doping concentration. 
Since the competing interactions involved in the SOS case are the double-exchange 
and superexchange, both of which are the nearest-neighbor interactions (thus both 
are strong) , the phase transition temperature is expected to be higher than the spi- 
ral spin orders in undoped i?Mn03 which are stabilized by the weak next-nearest- 
neighbor interactions. According to the Monte Carlo simulation, the multiferroic 
Tq may reach 100 KP^The ferroelectric polarization of SOS phase has dual origins: 
the Dzyaloshinskii-Moriya interaction effect at the stripe boundaries and the ex- 
change striction effect within every stripes. Thus, the expected polarization should 
be considerable. 

Recently, Liang et al. generalized the SOS phase to other doping concentrations 
x = 1/N (N = 2, 3, 4, ... , oo)p3A general SOS^ family has been revealed, which 
corresponds the CxEi— x family. The multiferroic phase predicted by Giovannctti 
et al. can be included into this family as a special case SOS 1/2 P Some of these 
SOSa; and C^Ei^ are multiferroic. More interestingly, there are various derivative 
superstructures of SOS^ and C^Ei-^, which are energy-degenerated with the cor- 
responding primary ones. In other words, there is dimensional reduction in these 
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Fig. 14. Spin patterns of the (a) SOS and (b) C-i/^E^^ phases, (c) The (common) spin structure 
factor of these two phases. The characteristic peaks appear at (1/4, 1/4), (1/8, 3/8), and (3/8, 
1/8) and their equivalent positions in the q + (1/2, 1/2) region. The ideal amplitudes of these 
three peaks are in the ratio of 2 : 1 : 1. Density-of-states of the (d) SOS and (e) C1/4E3/4 phases 
(at Jahn- Teller coefficient A = 0, fixing the spin pattern). Both of these phases have an energy gap 
at the Fermi level, implying the insulating fact. Reprinted figure with permission granted from S. 
|Dong et al, Phys. Rev. Lett. 103 (2009) 1Q7204[ M1 Copyright ©(2009) by the American Physical 
Society. 



complex spin structures.^ 



3.2. CaMn 7 12 

Although there are many theoretical predictions for new multiferroics, few of them 
have been verified experimentally. Some predicted materials are very difficult to syn- 
thesize, or need very rigorous conditions. For example, to obtain the aforementioned 
SOS phase, one needs to get a quarter-doped perovskite manganites with very nar- 
row e g bandwidth but weak Jahn- Teller distortions.^ In normal perovskite man- 
ganites, the bandwidth decreases with decreasing size of A-site cations. However, 
the Jahn- Teller distortions will be enhanced simultaneously. Therefore, it seems 
almost impossible to simultaneously fulfill both conditions in normal perovskite 
manganites. Furthermore, the quenching disorder of A-site cations will seriously 
suppress the fragile long-range SOS order. This disorder problem does not exist in 
undoped i?Mn03. However, in doped manganites, the size and valence differences 
between the rare-earth cations and alkaline-earth cations are crucial to determine 
physical properties F^3 Sr 2+ is certainly too large to reduce the bandwidth, while 
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Fig. 15. Crystal structures of quadruple perovskite CaMn70i2- Left: cubic Im3 at high tempera- 
tures (> 440 K). Right: rhombohedral i?3 at low temperatures (< 440 K). Green Mn's are at the 
A-site. Drawn by VESTA.ED 

there is no size- matched rare-earth cation for Ca 2+ . 

CaMn70i2, a quadruple perovskite, provides a unique structure to realize com- 
plex spin structures, e.g. the SOS phase. As shown in Fig. [l5j its lattice is cubic at 
high temperatures, with four perovskite units in each cell. This cubic cell distorts 
into a rhombohedral lattice at ~ 440 K, accompanying a charge-order transition.^ 
In CaMn70i2, Ca 2+ and three Mn 3+ occupy the A-site, while the rest four Mn's 
at the B-site have an average valence +3.25, which just corresponds to the quarter- 
doping. Comparing with normal perovskite manganites, CaMn70i2's lattice size is 
greatly shrunk and the B-site Mn-O-Mn bonds are seriously distorted due to the 
small A-site Ca 2+ and Mn 3+ , as required to reduce the e g bandwidth. Furthermore, 
the Jahn- Teller distortions are very weak in CaM^Oi^J^as also required to stabi- 
lize the SOS phase. The last but not the least, the A-site cations are fully ordered 
in quadruple perovskites, which eliminates the quenching disorder in normal doped 
perovskites. However, in the past years, CaMn70i2 have not been much studied 
and no direct evidence of multiferroicity has been reported. 

In 2011, Zhang et al. synthesized polycrystalline samples of CaMn70i2, and 
measured the ferroelectric polarization and magnetoelectric responses!^ As shown 
in Fig. [16} the ferroelectricity is quite remarkable. First, the ferroelectric transition 
temperature is much higher comparing with other magnetic multiferroic mangan- 
ites: T c = 90 K, which coincides with the Neel transition T N1 . Second, the po- 
larization reaches ~ 450 /iC/m 2 at 8 K when using a poling field 7 kV/cm, which 
seems to be unsaturated yet. Third, this polarization can be significantly suppressed 
by ~ 30% under a 9 T magnetic field, which also suggests the magnetic origin of 
polarization. 

Later, Jahnson et al. synthesized single crystal samples of CaMn.7012 and deter- 
mine the magnetic structures using the neutron powder diffraction.^ The polariza- 
tion, along the c-axis of the rhombohedral cell (the [111] qseudocubic axis), reaches 
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Fig. 16. (a) The pyroelectric currents of CaMri70i2 with a poling field E = 7 kV/cm. An anoma- 
lous contribution appears around Tn2 = 48 K (another magnetic transition temperature which's 
detail is unclear). Inset: The corresponding ferroelectric polarization around Tn2- (b) The ferro- 
electric polarization at 8 K as a function of the poling field E, which does not saturate up to 7 
kV/cm. (c) The suppression of polarization P(H) by a magnetic field H. Inset: the magnetoelectric 
response ratio, defined as [P(0) — P(H)]/ P(0) X 100%. Reprinted figure with permission granted 
from |G. Zhang et al, Phys. .Rev. B 84 (2011) 174413^ 1 Copyright ©(2011) by the American 
Physical Society. 



2870 /iC/m 2 at low temperatures, which is one of the largest measured polarizations 
in magnetic multiferroics, as shown in Fig. |17[ a). Most importantly, the magnetic 
structure has been resolved. As shown in Fig. |17[ c), the magnetic wave vector is (0, 
1, k), where k is ~ 0.963 between 90 K and 48 K but splits into two branches below 
48 K (~ 0.88 and ~ 1.042 at 5 K). The in-plane and out-of plane spin patterns are 
shown in Fig. 17 d,e). The incommensurate wave number k suggests a screw type 
spiral propagating along the c-axis, which is different from the cycloidal spiral in 
undoped -RM11O3. Within the ab plane, spins are also noncollinear. 

However, the complex spin order in CaMn 7 0i 2 is not the expected SOS phase, 
which may be due to the magnetic interaction between the A-site Mn and B-site 
Mn. In the aforementioned model calculation, the A-site magnetism have not been 
taken into consideration, which needs further theoretical investigations. 

Then, how to understand the giant improper ferroelectricity in CaMn70i2? 
According to the formula of Dzyaloshinskii-Moriya interaction (or the KNB theory), 
the screw type spiral along the c-axis can not give rise to a polarization, as in the Y- 
type hexaferrite Bao.sSri.sZ^Fe^C^-'^' Johnson et al. proposed that the ferroaxial 
coupling mechanism is responsible to the polarization. The symmetry group, i?3, 
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Fig. 17. (a) Polarization of CaMii70i2 single crystal along the hexagonal c axis (black) and the 
three qseudocubic < 100 > axes. Inset: A small anomalous at Tnj2- (b) The magnetic susceptibility 
parallel and perpendicular to the hexagonal c axis, measured under a 500 Oe magnetic field, in both 
zero field cooled and field cooled (500 Oe) conditions. Inset: the in-plane magnetic susceptibility 
at Tni. (c) Temperature dependence of incommensurate magnetic propagation along the c-axis. 
(d) and (e) are the magnetic structure of CaMn70i2 in the ab and ac planes, respectively. Mnl, 
Mn2, and Mn3 are shown in red, black, and yellow, respectively. The moments rotate in the ab 
plane with a circular envelope as depicted in (e). Reprinted figure with permission granted from 
|R. D. Johnson et al, Phys. Rev. Lett. 108 (2012) 06720l[ MI Copyright ©(2012) by the American 
Physical Society. 



belongs to the 3 ferroaxial class, for which a homogeneous structural rotation exists 
with respect to the high-temperature cubic group, represented by an axial vector 
A. Symmetry considerations allow a possible coupling between the axial vector A, 
the spin chirality a (= ■ (Si X Sj), which is a time-even, parity-odd pseudoscalar) 
along the c-axis, and the polarization P c along the c-axis. This P c ctA term in free 
energy will give rise to a finite P z when a is finiteEEH 

However, a recent density functional theory calculation showed that the screw 
spiral along c-axis is not crucial to obtain the giant ferroelectric polarization in 
CaMn70i2. ? Even with a zero c, the in-plane noncollinear spin pattern can induce 
a polarization up to ~ 4500 /iC/m 2 along the c axis. The underlying mechanism 
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are the combination of Dzyaloshinskii-Moriya interaction and exchange strictionP^ 
Therefore, further theoretical and experimental studies are needed to clarify the 
multiferroic mechanisms involved in CaMn70i2- 

4. Summary and perspective 

In this Brief Review, we have introduced some recent theoretical and experimental 
progress on multiferroic perovskite manganites, in which the ferroelectric polariza- 
tions are induced by particular magnetic orders. This review have (partially) em- 
phasized three fundamental questions: 1) Why and how can the spiral spin order 
and E-type antiferromagnetic order generate improper ferroelectric polarizations? 
2) Why are there the spiral spin order and E-type antiferromagnetic order in per- 
ovskite manganites and how to tune the competition between them? 3) Are there 
new magnetic multiferroics in doped manganites with better performance and how 
to understand their exotic properties? 

To answer these questions, we have taken the undoped i?MnC>3 as the example. 
The Dzyaloshinskii-Moriya scenario (or the Katsura-Nagaosa-Balatsky theory) and 
exchange striction effect have been applied to explain the origin of ferroelectricity. 
In some cases, these two mechanisms coexist and are mutually coupled. Both the 
classical spin model and quantum double-exchange model have been introduced 
to understand the phase diagram of i?Mn03. The theoretical SOS phase and ex- 
perimental CaMn70i2 have also been introduced as new magnetic multiferroics in 
doped manganites. 

The research field of multiferroicity remains very active and are developing very 
fast. The predictions/discoveries of new multiferroics have never been stopped since 
2003. Accompanying with the progress on new multiferroic materials and improved 
samples' quality, the multiferroic performance has been improved gradually, includ- 
ing both the working temperatures and ferroelectric polarizations. In addition, the 
theoretical understanding of multiferroicity has also been pushed forward in the 
past years. Even though, there remains many puzzles and debates, some of which 
have also been presented in this review. Hope future experimental and theoretical 
works will solve them. 
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